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■ Abstract. We have investigated the baryon-mass content in a subsample of 19 clusters of galaxies extracted from the X-ray 

flux-limited sample HIFLUGCS according to their positions in the sky. For these clusters, we measured total masses and 
characteristic radii on the basis of a rich optical spectroscopic data set, the physical properties of the intracluster medium (ICM) 
I . using XMM-Newton and ROSAT X-ray data, and total (galaxy) stellar masses utilizing the SDSS DR7 multi-band imaging. The 

O ' observed (hot) gas-mass fractions are almost constant in this mass range. We confirm that the stellar mass fraction decreases as 

, the total mass increases and shows (20 ±4)% scatter; in addition, we show that it decreases as the central entropy increases. The 

^ ' latter behavior supports a twofold interpretation where heating from merging quenches the star-formation activity of galaxies in 

massive systems, and feedback from supemovae and/or radio galaxies drives a significant amount of gas to the regions beyond 
rsoo or, alternatively, a substantially large amount of intracluster light (ICL) is associated with galaxies in nonrelaxed systems. 
04 ' Furthermore, less massive clusters are confirmed to host less gas per unit total mass; however, they exhibit higher mass fractions 

^ \ in metals, so that their ICM is more metal-rich. This again supports the interpretation that in the potential wells of low-mass 

systems the star-formation efficiency of galaxies was high or, alternatively, some gas is missing from the hot phase of the ICM. 
The former hypothesis is preferred as the main driver of the mass-dependent metal enrichment since the total mass-to-optical 
luminosity ratio increases as the total mass increases. 



o 

cn 

o 



On ' Key words. Cosmology: observations — Galaxies: clusters: general — Methods: data analysis — Surveys — X-rays: galaxies: 

»' . clusters — Galaxies: stellar content 



- ' — \ 

X 



1. Introduction /W»/Mtot, the stellar-to-total mass ratio in which the stars in 

galaxies are considered). We can use X-ray data alone to mea- 

Galaxy clusters are the largest gravitationally bound structures ^^^^ g^^.^^^^ fraction (e.g., Ettori et al. 2002; Zhang et 

^ , in the Universe, and they have become an important cosmo- 2006). There are two well-established methods of measur- 

- - ■ logical probe to constrain dark energy (e.g., Vikhlinin et al. .^^^ ^^^^j ^^^jj^^ ^^^^ ^ ^^j^^^ ^j^^^^^ 

2009). Calibrating the baryon content of galaxy clusters is thus »u » n * »u » » 1 1 w * 1 

^ & .7 b J g^j-g j-jjg stellar mass tor the total galaxy population ot a clus- 

a key input for understanding the structure-formation history, . f . » 1 1 • u j »u 1 1 • 

^ ° -" ter from the total lurmnosity based on the galaxy luminosity 

as well as systematic uncertainties in constraining cosmology ^^^^^-^^ (qlF) and a stellar mass-to-optical light ratio that de- 

with galaxy clusters as standard candles. p^^^^^ ^^^^ ^ ^j^^^^^ ^^^.^ j^^q. ^^^^^^j 

The hot intracluster medium (ICM; e.g., Sarazin & Bahcall 200O; Lin et al. 2003). The other is to measure the stel- 

1977) and the stars bound to member galaxies have been con- ^^^^^^ ^^e individual member galaxies and to construct 

sidered to be the main contributors to the baryon content in ^^e stellar mass function in order to sum these stellar masses 

clusters. The baryon-mass fraction of a galaxy cluster was thus ^^^^-^^ 3 q-^^-^- 2009). The contribution to 

defined as the sum of the gas-mass fraction (/g^s := M.^,JM,,„ ^^e total stellar mass budget associated with intracluster light 

the gas-to-total mass ratio) and the stellar mass fraction (/, ;= (j^L; e.g., Zibetti et a. 2005) could be significant especially 

on the mass scale of galaxy groups if the ICL makes more 

Send offprint requests to: Y.-Y. Zhang ^^^^ ^q_2q^^ ^^^^j U ^ Gonzalez et al. 2007). 
** Guest astronomer at the MPE 
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The discrepancies between measurements of the baryon-mass 
fraction of galaxy clusters and groups (e.g., Lin et al. 2003; 
Gonzalez et al. 2007; Krick & Bernstein 2007; Lagana et al. 
2008, 2011; Giodini et al. 2009; Andreon 2010; Dai et al. 
2010) and the measurement from the Wilkinson Microwave 
Anisotropy Probe (WMAP) 5-year data (Dunkley et al. 2009) 
allow for an ICL-to-total stellar mass ratio in a wide range up 
to ~ 50% with a controversial mass dependence. Most observa- 
tional studies do not support dominant contributions from the 
ICL (e.g., Zibetti et al. 2005; see Arnaboldi & Gerhard 2010 for 
a recent review). Furthermore, theoretical studies reach contro- 
versial conclusions on the origin and mass dependence of the 
ICL fraction (e.g., Murante et al. 2004; Dolag et al. 2010 and 
references therein). 

Understanding the baryon-mass fraction and its mass de- 
pendence will be a milestone to understand astrophysics in 
galaxy clusters, e.g., the origin of the ICL (e.g., Pierini et 
al. 2008), star-formation history (e.g., Fritz et al. 2010), and 
metal-enrichment history (e.g., Kapferer et al. 2009). This also 
helps to control systematic uncertainties in high-precision clus- 
ter cosmology. 

The HIFLUGCS (Reiprich 8l Bohringer 2002) provides an 
X-ray flux-limited sample of 64 nearby {z < 0.1) clusters se- 
lected from the ROSAT All-Sky Survey (RASS). Nineteen clus- 
ters in the HIFLUGCS have been imaged beyond 9r50(ll] in the 
Sloan Digital Sky Survey (SDSS, DR7, Abazajian et al. 2009). 
Since these 19 clusters were selected from the HIFLUGCS 
based on their sky positions, this subsample is unbiased with 
respect to the HIFLUGCS. Our sample has therefore inherited 
the nature of the flux selection. 

In this paper, we analyze XMM-Newton and ROSAT X-ray 
data, as well as optical data for this sample of those 19 clus- 
ters, and investigate the baryon content in clusters and its mass 
dependence. This study, similar to Andreon (2010), has the ad- 
vantage over many previous studies that the cluster total masses 
and reference radii are derived from the velocity dispersion 
based on optical spectroscopic data (see Zhang et al. 20 11; and 
references therein) that are independent of both the X-ray data 
used to determine the gas masses and the optical imaging data 
used to measure the stellar masses. It enables a relatively objec- 
tive investigation of the intrinsic scatter in the gas-mass, stellar 
mass and baryon-mass fractions including contributions from 
the intrinsic scatter of mass estimates. 

We organize the paper as follows. In Sect. |2] we describe 
the sample and data analysis. We present our results in Sect. [3] 
as well as discussions in Sect. |4l and summarize our find- 
ings in Sect. |5] We assume Q.^ = 0.3, Qa = 0.7, and Hq = 
70km s ' Mpc ' . Throughout this paper we apply the BCES re- 
gression fitting method taking measurement errors into account 
in both variables and their covariance (Akritas & Bershady 



' The cluster radius, r^, is the radius within which the mass den- 
sity is A times the critical density, pdz) = E^(z)3H^/(SnG), where 
E'^iz) = fiA + £!,„( 1 -H z) V ( 1 - - n,„)( 1 + zf . In this work, we mea- 
sured two cluster radii, r^m and r2m- As shown in Lagana et al. (2011), 
the coverage of 9r5oo is required to perform a robust background sub- 
traction to construct the GLF of the cluster. 



1996). Confidence intervals correspond to the 68% confidence 
level. 

2. Sample and data 

2.1. Sample 

The HIFLUGCS is an X-ray flux-limited sample of 64 nearby 
(z < 0.1) clusters selected from the RASS. Nineteen clusters 
(Table [B in the HIFLUGCS have available XMM-Newton and 
SDSS DR7 data (Abazajian et al. 2009). 

2.2. Cluster total mass and radius 

Recent observations suggest deviations from hydrostatic equi- 
librium in galaxy clusters (e.g., Zhang et al. 2008, 2010; 
Mahdavi et al. 2008). We therefore avoid using the X-ray hy- 
drostatic mass (e.g., Zhang et al. 2009) or caustic mass (e.g., 
Rines & Diaferio 2006) as the cluster mass. Furthermore, to 
guarantee a fair study, we choose to obtain the total mass from 
a quantity that is not linked to either quantities that we study, 
i.e., stellar mass and gas mass. 

We calculate the cluster mass from the velocity dispersion 
using the optical spectroscopic data listed in Table 1 in Zhang 
et al. (2011; also see references therein), which are independent 
of both the X-ray data used to derive the ICM properties and the 
optical SDSS imaging data used to derive the stellar masses. 

Taking the "harmonic" velocity dispersion (cr^ p) measured 
within an aperture of 1.2 Abell radii (a = 2.57 Mpc) from 
Zhang et al. (20 11), we follow the method described in Sect. 3 
in Biviano et al. (2006) of computing the cluster mass Ma-. 
We first determine an initial estimate of the mass from V3cra,p 
through eq. (2) in Biviano et al. (2006). An estimate of the ra- 
dius fv is derived through steps 7 and 9 in Biviano et al. (2006). 
Replacing the true quantities and cr^ with their estimates r\ 
and VScTg p in Fig. 4 in Biviano et al. (2006), we obtain an esti- 
mate of cTy. The cluster mass M^- is given by eq. (2) in Biviano 
et al. (2006) and is accurate to a; 10%. The mass error is derived 
by combining in quadrature the error in the velocity dispersion 
converted in mass, additional velocity dispersion error intro- 
duced by the uncertainty in Figure 4 in Biviano et al. (2006), 
and the 10% accuracy in Ma-- Typical mass errors are ~20% 
(Table [TJ. To compute r^oo from M200, we assume an NFW 
(Navarro et al. 1997) model with the concentration parameter 
given in step 7 in Biviano et al. (2006). 

As shown in Figure 3 in Zhang et al. (2011), the 
HIFLUGCS itself is biased toward X-ray luminous clusters at 
the low-mass end and contains an increasing fraction of merg- 
ers toward the high-mass end. This leads to a large fraction 
of X-ray luminous clusters in the low-mass regime and a domi- 
nant fraction of mergers in the high-mass regime for the sample 
of 19 clusters. The velocity dispersion estimates can be biased 
high by substructures, which leads to overestimating the dy- 
namical masses for mergers (e.g., Biviano et al. 2006). On the 
other hand, the cluster masses, derived from our best fit of the 
X-ray hydrostatic mass versus gas mass relation of the clusters 
and groups in Vikhlinin et al. (2006), Arnaud et al. (2007), Pratt 
et al. (2009), and Sun et al. (2009), are biased high for the X-ray 
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luminous clusters in the low-mass regime. Consequently, we 
observe that, at the low-mass/high-mass end, the total masses 
derived from the velocity dispersion appear low/high compared 
to the total masses derived from the gas masses as shown in 
Figure [T] The two most outliers are A2029 and A2065, which 
have the highest masses derived from the "harmonic" velocity 
dispersion. Therefore, the full sample of the 19 clusters and a 
sample of 17 clusters excluding A2029 and A2065 will be in- 
vestigated respectively, in Sect. [3] 

We derived all quantities consistently within the cluster ra- 
dius, rjoo, unless explicitly stated otherwise. 

2.3. X-ray data 

Details on the joint X-ray XMM-Newton and ROSAT data anal- 
ysis and properties of the XMM-Newton observations can be 
found in Zhang et al. (2009, 201 1). The procedures for clean- 
ing the data are described in Sect. 2.2 and for detecting and sub- 
tracting point-like sources in Sect. 2.3 in Zhang et al. (2009). 
The background treatment can be found in Sect. 2.4 of Zhang 
et al. (2009). Significant substructure features clearly detected 
in the image are excised before we perform the spectral and 
surface brightness analysis. 

The truncation radii of the XMM-Newton observed surface 
brightness profiles, within which the signal-to-noise ratio (S/N) 
is greater than 3, are rather small (< rsoo). The ROSAT ob- 
served surface brightness profiles cover radii well beyond rjoo 
with S/N > 3, although with sparse data points to resolve 
the cluster core. As noted in Sect. 2.2 in Zhang et al. (2011), 
we directly convert the ROSAT surface brightness profile to 
the XMM-Newton count rate using the best-fit spectral model 
obtained from the XMM-Newton data. We then combine the 
XMM-Newton surface brightness profile within the truncation 
radius, where the XMM-Newton S/N is ~ 3, with the ROSAT 
converted surface brightness profile beyond the truncation ra- 
dius for further analysis. 

The soft band (0.7-2 keV) X-ray surface brightness profile 
model SxiR), in which R is the projected radius, is linked to 
the ICM radial electron number density profile ne(r) and emis- 
sivity function as an integral performed along the line of sight 
convolved with the XMM-Newton point-spread function (PSF) 
as a function of the detector coordinates {X, Y), 

SxiR)'^ J J J PSF{X,Y,R)nl[^[F^^didXdY. (1) 

The gas mass distribution is derived by integrating the gas den- 
sity distribution. Gas masses are measured with respect to rjoo 
derived in Sect. 12.21 Errors are derived by combining in quadra- 
ture measurement errors of gas masses at rsoo and errors due to 
the variation of rsoo within its Icr interval. Typical gas mass 
errors are within 10% (see Table[TJ. 

The spectral analysis through which we obtain the radial 
temperature distribution T{r) is documented in Sect. 3 in Zhang 
et al. (2009). The radial entropy distribution is computed as 
S{r) - T(r)n^^^^(r). The XMM-Newton data enable measuring 
the radial temperature distribution beyond 0.2r5oo. The central 
entropies are thus measured at O.lrsoo and 0. 2^500. 



2.4. Optical imaging data 

The SDSSH DR7 has the advantage over previous SDSS data 
releases of using more reliable sky background subtraction 
(Abazajian et al. 2009). We used the magnitudes in the "dered" 
table in the "PHOTO" catalog. Hereafter we briefly summarize 
the procedure on constructing the GLF of the galaxy cluster 
(see Lagana et al. 2008, 201 1 for details) as follows. 

A color-magnitude diagram is constructed for the galax- 
ies within a projected radius of rsoo. The galaxies with 0.7 < 
(g -i) < 1 .6 and / < 1 8 in the g-ivs. i diagrarr|| are bright and 
red galaxies, and are used to make a preliminary fit of the red 
sequence in the color-magnitude diagram with a linear relation. 
Iteratively, the red sequence is fitted using those galaxies in the 
full magnitude range within 3cr redward and blueward of the 
previous red-sequence fit until the fit converges. The red galax- 
ies are defined as those within 0.3 mag redward and blueward 
m g - i of the best-fit red sequence, and blue galaxies as those 
bluer than the lower limit for the red galaxy selection. 

The fore- and background galaxies are subtracted statisti- 
cally (e.g., Zwicky 1957; Oemler 1973). The /-band GLF is 
constructed for the selected red and blue galaxies with a bin 
size of 0.5 mag. The GLF of the background galaxies is de- 
rived using those galaxies within an annular region located 
beyond 8r5()o from the X-ray flux-weighted cluster center A 
power-law model is fitted to the GLF of the background galax- 
ies and subtracted from the overall GLF. A fair background 
subtraction should not introduce major uncertainties since the 
GLF is sampled well enough at the bright end, and it does not 
abruptly rise at the faint end. Large-scale structures (LSSs) may 
introduce some uncertainties into this kind of background sub- 
traction. Recent tests (e.g., Paolillo et al. 2001; Gonzalez et 
al. 2007, Filippis et al. 2011) show that the difference in the 
GLF is within 10% compared to the case using the background 
level derived from the SDSS GLF (e.g., Blanton et al. 2003), 
integrating over redshift and applying mean evolutionary and 
/T-corrections (e.g., Bruzual & Chariot 2003) with the Padova 
isochrones (Bertelli et al. 1994). Corrections due to poor or ab- 
sent sampling of the faint end of the GLF should only be a few 
per cent (see Andreon 2010). 

There are two obvious components in the GLF for most 
clusters. A double Schechter (1976) function is thus used to 
fit the GLF of each cluster, in which the covariance between 
the power-law index at the faint end {a) and the characteris- 
tic magnitude (M*) is taken into account (see Table |2]i. The 
/iT-correction for early-type galaxies and the redshift evolution 
correction given in Poggianti (1997) are adopted to correct the 
characteristic magnitudes. The total /-band luminosity is de- 
rived by integrating the double Schechter function for magni- 
tudes brighter than -14. 

^ The A'-band luminosity from the Two Micron All Sky Survey 
(2MASS) observations can also be used to measure the stellar mass 
(Lin et al. 2003). Compared to the 2MASS observations, the SDSS 
data allow us to make a more reliable selection of cluster member 
galaxies, which greatly reduces fore- and background contaminations. 

' This magnitude cut best suits the majority of the 19 clusters. For 
four clusters (IIIZw54, A1367, MKW4, and A400) we apply / < 16, 
which better defines the preliminary slope of the red sequence. 
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The mass-to-light ratios for ellipticals and spirals, respec- 
tively, following Kauffmann et al. (2003) assuming the Salpeter 
(1955) initial mass function (IMF) are adopted to compute the 
stellar mass from the optical luminosity. The stellar mass es- 
timate is tied to the choice of IMF. Changing the IMF will 
scale the stellar mass estimate by a fixed factor, e.g., from a 
Kroupa (2001) IMF to a Salpeter (1955) IMF with a cut-off at 
O.IMq results in a factor of 2 increase in the stellar mass (e.g., 
Kauffmann et al. 2003). 

Any systematic offset in the spectrophotometric calibration 
will change the stellar mass estimate (e.g., Kauffmann et al. 
2003). Fritz et al. (2010) have carried out a thorough compar- 
ison of stellar mass estimates for a large sample of more than 
50 nearby clusters (z < 0.07) computed from spectroscopic 
and photometric data. They found that the stellar masses com- 
puted from the fiber-aperture magnitude using SDSS DR7 pho- 
tometric data are lower by ~ 0. 15 dex than the values computed 
from the spectroscopic data within the same aperture (also see 
Andreon 2010). Our sample is in similar redshift and mass 
ranges as their sample, and we also calculated the total stel- 
lar mass using the SDSS DRV data. Therefore, this bias is cor- 
rected in our total stellar mass estimates. 

3. Results 

Figure |2] shows the gas-mass, stellar mass, and baryon-mass 
fractions as a function of the total mass. We also compared our 
observational sample to some existing observational samples 
(e.g., Lin et al. 2003; Gonzalez et al. 2007; Giodini et al. 2009; 
Pratt et al. 2009; Sun et al. 2009; Andreon 2010; Dai et al. 
2010), as well as simulated samples (e.g., Puchwein et al. 2008, 
2010; Fabjan et al. 2010). 

3.1. Gas-mass fraction 

The gas-mass fractions of the 19 clusters appear flat in the 
given narrow mass range with (26 + 4)% intrinsic scatter. 
Excluding the two most prominent outliers, A2029 and A2065, 
shown in Sect. 12.2! the gas-mass fraction increases with the to- 
tal mass as /gas.500 = lO-" '"^°"''(M5oo/[lO''*M0])"-3'^±°-3'' with 
(23 + 17)% intrinsic scatter To constrain the slope precisely, 
we combine our 19 clusters and the groups in Sun et al. (2009), 
which show an increasing gas-mass fraction as a function of the 
total mass, i.e., /gas.soo = lO-(' °^=^°°2)(M5oo/[lO'4M0])°-'"'=too^ 
with (26 ± 8)% intrinsic scatter. This shows the importance of 
having a wide mass range to calibrate the mass dependence of 
the gas-mass fraction. 

The normalization for our sample is comparable to what is 
obtained for the X-ray selected cluster sample in Pratt et al. 
(2009) and X-ray selected group sample in Sun et al. (2009), 
and is higher than the 2MASS selected stacked groups in Dai 
et al. (2010) with > 3cr significance. This indicates that to 
some extent X-ray selected groups may be biased toward sys- 
tems with high gas-mass fractions, whereas near-infrared (or 
optically) selected groups are biased toward systems with low 
gas-mass fractions. (The extreme systems are probably the 
so-called X-ray under-luminous systems, e.g., Dietrich et al. 
2009.) 



The gas-mass fractions of our observational sample agree 
well with the predictions of the simulated sample, in particular, 
in the high-mass regime, i.e., > 1O''^M0. It is interesting that the 
slope of the mock sample with AGN feedback is between the 
measured slopes of X-ray selected and near-infrared selected 
samples. This can either be because that the simulation sample 
is free of selection biases, or our sample has a sampling bias in 
the low-mass regime caused by low-quality statistics. 

3.2. Stellar mass fraction 

The stellar mass fraction in our observational sample 
decreases with increasing cluster mass; i.e., /*,5oo - 
^0-(i.53±o.05)(^^^^/^Qi4]y^^])-(0.49±o.09) ^jjjj + 4)% intrin- 
sic scatter for the 19 clusters. Excluding A2029 and A2065, the 
best fit is /.,5oo = 1O-(1-5«±oo5>(M5oo/[1O'4M0])-(O-3''±"'») with 
(15 ± 29)% intrinsic scatter We also derived the local regres- 
sion fit used in Zhang et al. (2009), which demonstrates the 
flattening of the stellar mass fraction in the mass range lower 
than 2x IO'^^Mq (see Fig.|2]i. There are three situations that may 
introduce measurement uncertainties and affect the best fit, (1) 
the poor sampling in the low-mass regime, i.e., < lO'^M©, (2) 
the potential volume dependence in the background estimate, 
and (3) the two outliers for which the total masses are likely 
overestimated. The volume fraction of the background galax- 
ies decreases with decreasing cluster redshift. The estimate of 
the local background in deriving the GLF becomes less reliable 
in, e.g.. Coma (rsoo = 41.1'). 

Within Icr, the best fit of the stellar mass fraction as a 
function of the total mass agrees between our sample and the 
published samples (e.g., Lin et al. 2003; Gonzalez et al. 2007; 
Giodini et al. 2009; Andreon 2010; Dai et al. 2010) as shown 
in Table [3] The stellar mass fractions of nearly half of the 19 
clusters in the observational sample agree with those of the 
simulated samples within their Icr interval (e.g., Fabjan et al. 
2010; Puchwein et al. 2010). The amount of the stellar mass 
in galaxies in simulations depends rather sensitively on how 
the supernova feedback is included (e.g., Springel & Hernquist 
2003a; Borgani et al. 2006). For example, using kinetic feed- 
back, e.g., in addition to the thermal feedback that Puchwein et 
al. (2010) used, one can obtain up to a factor of 2 fewer stars. 
Furthermore, the cluster mass is well defined in simulations, 
while there are known uncertainties in mass measurements of 
galaxy clusters (e.g., Biviano et al. 2006). This also causes dif- 
ferent trends and scatter between the observational and simu- 
lated samples. 

3.3. Baryon-mass fraction 

In simulations, part of the gas does not enter the potential well 
of a galaxy group because gas is collisional. This phenomenon 
drives the so-called gas depletion. This "correction" is usu- 
ally taken into account when estimating the baryon-mass frac- 
tions for observed groups. It amounts to roughly 10% within 
rsoo for groups (e.g., Frenk et al. 1999; Kay et al. 2004). 
The gas depletion is only a few percent in galaxy clusters as 
massive as those considered in this study (e.g., Ettori et al. 
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2006; Fabjan et al. 2010). We therefore do not account for 
it, and define the baryon-mass fraction of the galaxy cluster 
as the sum of the stellar mass fraction and gas-mass fraction, 
/b,500 := M..500/M500 + Mgas.soo/MscK). 

In Fig. |2] we also show the baryon-mass fraction as a 
function of the total mass. The baryon-mass fractions of all 
19 clusters appear flat as a function of the total mass. Even 
when A2029 and A2065 are excluded, the best fit, /b.soo = 
lQ-(0.964±o.250)(^^^^j/^Oi'*Mo])°-224±0-567^ Still appears flat as a 

function of the total mass within its Icr error. 

As shown in Table|3] the baryon-mass fractions of our sam- 
ple display a similar trend as that exhibited by other obser- 
vational samples (e.g., Lin et al. 2003; Gonzalez et al. 2007; 
Giodini et al. 2009; Andreon 2010). According to the best fit 
for our sample, the value reaches the WMAP 5-year measure- 
ment of the cosmic baryon-mass fraction (Dunkley et al. 2009) 
at ~ 7.7 X IO'^^Mq. The observational sample displays large 
scatter suffering from measurement uncertainties and probably 
also cluster physics. 

4. Discussions 

4.1. ICL 

The small discrepancy between the baryon-mass fraction of our 
observational sample and the cosmic value appears to increase 
toward the low-mass end. The ICL is suggested to be one of 
the most important forms of missing baryons that accounts for 
this discrepancy, and accounts for 6%-22% of the total cluster 
light in the r-band (e.g., Krick & Bernstein 2007; Gonzalez et 
al. 2007; Zibetti et al. 2005). Pierini et al. (2008) suggest that 
the ICL fraction may increase in cluster mergers. 

Simulations with different cluster physics (e.g., Frenk et al. 
1999; Borgani et al. 2004; Kay et al. 2004; Kravtsov et al. 2005; 
Nagai et al. 2007; Evrard et al. 2008; Sijacki et al. 2007; Fabjan 
et al. 2010; Kapferer et al. 2010; Puchwein et al. 2010) give dif- 
ferent predictions. Kravtsov et al. (2005) performed numerical 
simulations that assume nonradiative hydrodynamics (without 
dissipations) in one run and radiative cooling, as well as sev- 
eral other physical processes in the other run. They find that 
the baryon-mass fraction does not depend on the cluster mass 
in nonradiative simulations, and that the predicted value within 
rsoo is 5% below the cosmic value. In simulations with radiative 
cooling, star formation, and AGN feedback, both the baryon- 
mass fraction and the gas-mass fraction increase with increas- 
ing mass (e.g., Kravtsov et al. 2005; Puchwein et al. 2010). In 
the runs with AGN feedback, the baryon-mass fraction is close 
to the cosmic value for the most massive clusters, but is sig- 
nificantly lower than the cosmic value for low-mass systems 
(e.g, Fabjan et al. 2010; Puchwein et al. 2010). This is due to 
the amount of the ICL, as well as of a significant amount of 
gas removed by AGN heating from the central regions of clus- 
ters and driven to > r5oo, in particular, in low-mass systems as 
shown by the sum of the baryon-mass fraction (/gas,500 + /«,50o) 
and the ICL in Figure 5 in Puchwein et al. (2010). 

In Puchwein et al. (2010), the ICL amounts to ~ 58% of the 
total stars, i.e., the sum of the ICL and stars in cluster galaxies, 
in a cluster of 4 x lO'^'M©. The amount of the stellar mass in 



the ICL found in simulations is enough to explain the difference 
between the baryon-mass fraction of our observational sample 
and the WMAP 5-year predicted cosmic value in the low-mass 
regime. However, such a large amount of the stellar mass in the 
ICL seems to be incompatible with recent observations (e.g., 
Krick & Bernstein 2007; Gonzalez et al. 2007; Zibetti et al. 
2005). This suggests that the ICL only accounts for part of the 
missing baryons and that other mechanisms like gas expulsion 
by AGN heating may be important. 

4.2. Star-formation efficiency 

We observe a decreasing stellar mass fraction as a function 
of the total mass (Fig. |2]i. The stellar mass fraction has been 
widely used as an estimate of the star-formation efficienc}0 
(e.g., Bryan 2000). There are no correlations between the 
stellar-to-gas mass ratio within rsoo and redshift. This indicates 
that the mass dependence of the star-formation efficiency is not 
a selection effect of redshift evolution. 

Other observational samples also show varying star- 
formation efficiency with total mass (e.g., Gonzalez et al. 2000; 
Linet al. 2003; Laganaet al. 2008, 201 1; Ettori et al. 2009). Dai 
et al. (2010) suggest that the fraction of baryon loss is deter- 
mined by the depth of the potential well of the system, which 
leads to an increasing baryon-mass fraction as a function of 
the increasing total mass of a cluster Toward the high-mass 
end, it approaches the cosmic value. Our finding supports two 
interpretations demonstrated in simulations (e.g., Fabjan et al. 
2010). In massive clusters, more hot gas and dark matter is set- 
tled in the deep potential wells than individual galaxies. In low- 
mass systems, the accretion of individual less massive galax- 
ies is more important, and more low entropy gas is brought in 
to form stars. Both processes result in low star-formation effi- 
ciency in massive systems. 

4.2.1 . Mergers and energy inputs from feedback 

In the hierarchical structure-formation scenario, the efficiency 
of stripping gas from member galaxies depends on the total 
mass of the system. In more massive clusters, galaxies exhibit 
higher velocity dispersion due to their deep potential wells, and 
as a consequence the efficiency of stripping gas from mem- 
ber galaxies is higher, which diminishes the star-formation 
efficiency in the galaxies. In low-mass systems, heating due 
to merging is less efficient because of their shallow potential 
wells. On the other hand, the effect of AGN feedback is signif- 
icant in groups due to their shallow potential wells. Giodini et 
al. (20 1 0) find that feedback from radio-mode AGN heating can 
account for the reduced gas fractions in groups in the COSMOS 
field. Simulations indeed show that a significant amount of gas 
has been removed by AGN heating from the central regions 
and been driven to the region beyond rsoo (e.g., Fabjan et al. 
2010; Puchwein et al. 2010). Furthermore, the fraction of X- 
ray luminous clusters is high in the low-mass regime because 

■* The star-formation efficiency used in this work is not the stellar- 
to-initial molecular hydrogen mass ratio, but the stellar-to-total cluster 
mass ratio. 
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of the nature of the flux selection of the HIFLUGCS. Most X- 
ray luminous systems are cool-core clusters and they host AGN 
in the central regions. There may be correlations between the 
stellar-to-gas mass ratio and the indicators of merging and var- 
ious energy inputs from feedback. 

As shown in Fig. [3] the scatter in the stellar-to-gas mass ra- 
tios appears to increase with decreasing off'set between the X- 
ray flux-weighted center and the bright cluster galaxy (BCG). 
We divided the sample into two comparable-size subsamples 
with a threshold of O.Olrsoo off'set, and found (41 + 8)% and 
(66 ±11)% intrinsic scatter for the ten high- and nine low- 
offset systems, respectively. This indicates that the scatter of 
the stellar-to-gas mass ratios may reflect the freedom in the 
formation history of the systems since the last major merging 
episode. 

Neither the stellar-to-gas mass ratio nor the offset between 
the X-ray flux-weighted center and the BCG correlates with 
the /-band magnitude difference between the BCG and the sec- 
ond brightest cluster galaxy. There are several possible expla- 
nations. The epoch of major star formation in the progenitors of 
massive ellipticals may happen much earlier than they fell into 
the cluster potential well. There may be smoother accretion of 
dark matter and hot gas than individual galaxies in quiescent 
systems. There may be a large fraction of the stellar mass in 
the ICL since mergers and tidal stresses have already happened 
and dynamical friction had more time to operate before the time 
of observations. 

There is no correlation between the stellar-to-gas mass ratio 
and the BCG magnitude. It indicates that the bulk of star for- 
mation does not dominate in the central region of the cluster, 
where the BCG is (also see Edwards et al. 2007; Loubser et al. 
2009; Wang et al. 2010). The star-formation efficiency may be 
similar in the cluster core and the outskirts. 

From the X-ray data we derived the cooling radius, rcooi, 
within which the cooling time is shorter than the age of the 
Universe. There is no evidence of high star-formation effi- 
ciency with increasing X-ray luminosity within the cooling 
radius. This indicates that radiative cooling losses within the 
cooling radii are offset by some heating mechanism(s). 

The stellar-to-gas mass ratio decreases with increasing 
entropy at O.lrsoo and 0.2r5oo (Fig.©, i.e., M,_5oo/Mgas,5()o = 

10-(0.477±0.075)5^^^^^^^/^00]^gV cm2])-(103rf-20) ^^^^ 

(51 + 7)% intrinsic scatter and Mt,5oo/Mgas,5oo = 
2Q-(0.i57±o.ii3)(5^^^^^^/^OOkeV cm^])-(i iv±o'i9) ^ith (41+5)% 
intrinsic scatter For the cool-core clusters alone, the scat- 
ter of the stellar-to-gas mass ratios becomes smaller, e.g., 
M,5oo/Mg,,,5oo = 10(2 i«-"-4«)(5o.i.3„o/[100keV cm']yii-'^-OA,, 
with (41 ±5)% intrinsic scatter The data support the interpre- 
tation that heating from merging quenches the star-formation 
activity of galaxies in massive systems, and feedback from 
supernovae and/or radio galaxies drives a significant amount of 
gas to the regions beyond rsoo or, alternatively, a substantially 
higher stellar mass fraction in the ICL is present in nonrelaxed 
systems (e.g., Pierini et al. 2008). 

Given the evidence of AGN feedback with radio-jet pow- 
ered energy input in nearby clusters (e.g., Birzan et al. 2008), 
we calculated the cavity power following the correlation in 
Birzan et al. (2008) from the radio bolometric luminosity given 



in Mittal et al. (2009). Unfortunately, we found no correlation 
between the stellar-to-gas mass ratio and the ratio between the 
X-ray bolometric luminosity within the cooling radius and the 
cavity power, respectively. The data do not support our guess 
on heating from the ongoing AGN feedback. 

The explanation may be that the cavity power provided by 
recent AGN activities only affects the current star-formation 
rate and it has less effect on more massive systems, in which en- 
ergy input from merging becomes more important. The stellar- 
to-gas mass ratio is set by the whole star-formation and energy- 
input histories of the galaxy clusters, and not by the current 
star-formation rates of the member galaxies. 

4.2.2. Metal enrichment 

The iron abundance of the ICM cannot be measured in a ho- 
mogeneous radial range for nearby clusters due to the small 
field-of-view (FOV) of XMM-Newton compared to the sizes 
of the clusters. This causes scatter in such measurements. On 
one hand, XMM-Newton FOV covers more of the cluster areas 
for lower-mass systems for a fixed redshift. Since the universal 
iron abundance decreases with increasing cluster radius (e.g., 
Kapferer et al. 2009, 2010), the measured global iron abun- 
dance within the XMM-Newton FOV is thus biased lower for 
low-mass systems than for high-mass systems. On the other 
hand, the X-ray emission of less massive systems is traced out 
to smaller fractions of rsoo than for more massive systems ac- 
cording to the flux limit. This should result in higher metallicity 
measurements in less massive systems when all systems can be 
observed out to large radii, e.g., ~ 0.5r5oo. We normalized the 
radius, within which the global temperature and metallicity are 
measured, by rsoo, and found that the normalized value has no 
correlation with rsoo. Therefore both effects are negligible for 
our sample. 

As shown in the left hand panel in Fig. |5j we 
found an increasing stellar-to-gas mass ratio with 
decreasing cluster mass, i.e., M, 5oo/M„as 500 - 
1O-*O-557±o<"I)(M5oo/[1O'4M0])-(O"7±o.ioi) ' ^ith °(29 + 5)% 

intrinsic scatter The cool-core (classified as "S" in Table [U 
clusters display a flatter slope, -(0.417 + 0.067), and smaller 
intrinsic scatter, < 17%, than the noncool-core (classified 
as "W" and "N" in Table [TJ clusters, which show a slope 
of -(0.617 ± 0.228) and (38 + 20)% scatter There is an 
increasing stellar-to-gas mass ratio with increasing iron 
abundance of the ICM in the right hand panel in Fig. |5j i.e., 

M,,5O0/Mgas,50O = 1O-("1O6±O.174)(2/Z0)1-53±O.38 ^jtjj (44 + 

intrinsic scatter We note that the trend holds even when one 
looks at the cool-core clusters and noncool-core clusters, 
respectively. It is interesting that the scatter of the cool-core 
clusters and noncool-core clusters is comparable. 

The above two correlations indicate that the iron in the ICM 
mainly comes from the pollution by the star formation that hap- 
pened in the past. In less massive galaxy systems, the star for- 
mation efficiency is higher. In other words, more stars were 
formed that have delivered more metals to enrich the hot gas. In 
high-mass clusters, energy feedback from, e.g., merging, may 
quench star formation in their member galaxies, which results 
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in low star-formation efficiency and less metal enrichment in 
the hot gas by stars. In addition, a larger amount of hot gas is 
accreted in a more massive system due to its deeper potential 
well, which dilutes the iron abundance more. 

As shown in Fig. |6] the observational sample shows de- 
creasing iron abundance with increasing total cluster mass, 
i.e., Z/Zo = 10-("-28^o.07)(^^^^/[1014Mo])-("-383±o.ioi) ^j^^ 

(38 ± 5)% intrinsic scatter, and increasing iron mass with 
increasing total cluster mass, i.e., ZMgas.soo/IZoMo] = 

jQl2.7±0.03(^^^^/^Ql4]y^^-|-)0.842±0.049 ^jtjj (17'+ 3)% intrinsic 

scatter. This finding agrees with the trend seen in simulations 
(e.g.. Figure 11 in Fabjan et al. 2010) that less massive clusters 
have lower gas-mass fractions but higher iron-mass fractions, 
so that the gas is more metal-rich in those systems. This can 
be explained in two ways. One is that the star-formation ef- 
ficiency is high in low-mass systems: more gas is converted 
into stars within rgoo which yields more metals that enrich 
the gas. Alternatively, metal-poor hot gas has not fallen into 
the potential wells in low-mass systems. Most likely, metal- 
loaded outflows from galaxies started well before the potential 
well of a cluster was settled (see Noll et al. 2009), so that the 
metal enrichment of the observed ICM reflects a heavier star- 
formation activity in low-mass systems. Accordingly, the to- 
tal mass-to-optical light ratio decreases with decreasing cluster 
mass (Fig.|7]l. 

4.2.3. Total mass-to-optical light ratio 

Optical and near-infrared observations suggest that the to- 
tal mass-to-optical light ratio, i.e., M/L™", increases with 
cluster mass with a slope in the range of 0.2-0.4 assum- 
ing a power-law relation (e.g., Adami et al. 1998; Bahcall 
& Comerford 2002; Rines et al. 2004). In Fig. H we show 
the total mass-to-optical light ratio as a function of the total 
mass for our sample, which follows a trend of Msoo/i^gyo = 

jq2.02±0.05(^^^j^/[1q14i^^-|-)0.468±0.085 ^jjj^ (^7 + 4)% intrinsic 

scatter. When A2029 and A2065 are excluded, the best fit is 

Msoo/Llf^ = 102"5=t005(^^^^/[jQl4]^^])0.371±0.098 ^jjjj ^3 + 

10)% intrinsic scatter 

The behavior of the total mass-to-optical light ratio sup- 
ports the variation in star-formation efficiencies found in sim- 
ulations (e.g., Springel & Hernquist 2003a; Saro et al. 2006). 
More massive clusters host more massive galaxies dominated 
by old, passively evolving stellar populations, where optical 
light has declined more significantly than the average co- 
eval galaxy due to their earlier formation time. Alternatively, 
more massive clusters may have accreted many more low-mass 
galaxies, or inhibited star formation in intermediate-mass/low- 
mass galaxies in a more efficient way, which leads to higher 
total mass-to-optical light ratios in massive systems. 

The scatter in the total mass-to-optical light ratio vs. total 
mass may indicate the coupling between the stochastic nature 
of the assembly times of dark matter halo and stellar mass. The 
best fit excluding A2029 and A2065 leads to Msoq/V;"^^ ~ 264 
for a system of 10 '^M©. Given the galaxy luminosity density 
and the assumption of an equivalent total mass-to-optical light 
ratio of the Universe to our extrapolated value for a M500 = 



IO'^Mq system, this results in a consistent Qn, with the concor- 
dance value within the scatter (also see e.g., Girardi et al. 2000; 
Reiprich 2003). 

5. Conclusions 

We have investigated the baryon content of a sample of 19 clus- 
ters selected from an X-ray flux-limited sample of 64 nearby 
clusters, in which we measured the cluster total masses and 
radii based on the optical spectroscopic data, which are inde- 
pendent of both the X-ray XMM-Newton and ROSAT data that 
we used to derive X-ray quantities, such as the gas masses, and 
of the SDSS DR7 imaging data that we used to compute the 
stellar masses. We summarize our results and interpretations as 
follows. 

- The gas-mass fraction of the observational sample in- 
creases with increasing cluster mass. The observed gas- 
mass fractions of our sample agree with previous observa- 
tional results (e.g.. Sun et al. 2007; Pratt et al. 2009). The 
values of our observational sample are also in good agree- 
ment with those of the simulated samples (e.g., Fabjan et 
al. 2010; Puchwein et al. 2010) in the high-mass regime. 

- The stellar mass fraction decreases with increasing clus- 
ter mass and exhibits large scatter (20 + 4%). The observed 
stellar mass fractions of our sample agree with previous ob- 
servational results and the results of the simulated clusters 
(e.g., Fabjan et al. 2010; Puchwein et al. 2010) within the 
scatter. 

- We observe lower stellar-to-gas mass ratios in those sys- 
tems that display higher central entropies indicating energy 
inputs from feedback. This supports a twofold interpreta- 
tion: heating from merging quenches the star formation in 
galaxies in massive systems, and feedback from supernovae 
and/or radio galaxies removes a significant amount of gas 
in low-mass systems or, alternatively, a substantially higher 
stellar mass fraction in the ICL is associated with dynami- 
cally active systems of galaxies. 

- Less massive clusters are confirmed to host less gas per 
unit total mass; however, they exhibit higher metal-mass 
fractions, so that their ICM is more metal-rich. Member 
galaxies of low-mass systems may be forming stars more 
efficiently, thereby producing more metals, part of which 
are ejected from the galaxies and enriched the ICM to a 
higher level. Alternatively, some gas could have been ex- 
pelled from these systems or be missing from the hot phase. 
The former hypothesis is preferred as the main driver of 
the mass-dependent metal enrichment since the total mass- 
to-optical luminosity ratio increases as the total mass in- 
creases. 
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Table 1. Properties of the 19 galaxy clusters 
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Name 




X-ray center (J2000) 
R.A. DEC 


redshift 


M500 
lO'+Mo 


10'"*] 


Mo 


10'^^ Mo 


M..500 

IO'^Mq 


Undisturbed 

/ ^^^^ 
1 Cool core 


A0085 


00: 


:41:50.306 


-09:18:11.11 


0.0556 


6.37 ± 


1.00 


5.68 ± 


0.37 


8.13 


± 0.38 


2.71 ± 1.83 


Y/S 


A0400 


02: 


:57:41.349 


-1-06:01:36.93 


0.0240 


1.83 ± 


0.39 


1.07 ± 


0.07 


1.36 


± 0.05 


1.58 ± 1.07 


N/N 


IIIZw54 


03: 


:41:18.729 


-1-15:24:13.91 


0.0311 


1.91 ± 


0.58 


1.18 ± 


0.08 


1.45 


+ 0.26 


2.79 ± 1.56 


Y/W 


A 1367 


11: 


: 44:44.501 


4-19:43:55.82 


0.0216 


1.76 ± 


0.27 


2.11 ± 


0.14 


2.07 


+ 0.07 


1.72 ± 0.30 


N/N 


MKW4 


12: 


: 04:27. 660 


-1-01:53:41.50 


0.0200 


0.50 ± 


0.14 


0.58 ± 


0.04 


0.47 


+ 0.02 


0.63 ± 0.18 


Y/S 


ZwC11215 


12: 


: 17:40.637 


-1-03:39:29.66 


0.0750 


4.93 ± 


0.98 


4.34 ± 


0.28 


6.10 


+ 0.29 


3.06 ± 2.30 


Y/N 


A1650 


12: 


:58:41.885 


-01:45:32.91 


0.0845 


3.44 ± 


0.66 


4.28 ± 


0.27 


5.09 


+ 0.73 


3.09 ± 1.43 


Y/W 


Coma 


12: 


:59:45.341 


-1-27:57:05.63 


0.0232 


6.55 ± 


0.79 


6.21 ± 


0.40 


8.42 


+ 0.63 


2.23 ± 1.12 


N/N 


A 1795 


13: 


:48:52.790 


+26:35:34.36 


0.0616 


3.41 ± 


0.63 


4.46 ± 


0.29 


5.11 


+ 0.14 


^ A A r\ a t 

2.44 ± 0.81 


Y/S 


MKWo 


14: 


: 40:42. 150 


+03:28:17.87 


0.0270 


0.62 ± 


0.12 


1.10 ± 


0.07 


0.80 


±0.12 


0.52 ± 0.39 


N/N 


A2029 


15: 


: 10:55.990 


+05:44:33.64 


0.0767 


14.70 d 


I 2.61 


6.82 ± 


0.44 


13.35 


±0.53 


4.49 ± 0.80 


Y/S 


A2052 


15: 


: 16:44.411 


+07:01:12.57 


0.0348 


1.39 ± 


0.28 


2.03 ± 


0.13 


1.86 


±0.10 


2.84 ±0.33 


Y/S 


MKW3S 


15: 


:21:50.277 


+07:42:11.77 


0.0450 


1.45 ± 


0.34 


2.29 ± 


0.15 


2.13 


±0.09 


1.30 ± 1.76 


Y/S 


A2065 


15: 


:22:29.082 


+27:43:14.39 


0.0721 


11.18d 


I 1.78 


3.35 + 


0.22 


7.66 


± 1.44 


3.04 ± 1.53 


N/W 


A2142 


15: 


:58;19.776 


+27:14:00.96 


0.0899 


7.36 ± 


1.25 


10.26 d 


:0.66 


13.76 


±0.73 


4.04 ± 2.22 


Y/W 


A2147 


16: 


:02: 16.305 


+ 15:58:18.46 


0.0351 


4.44 ± 


0.67 


3.63 ± 


0.23 


5.04 


±0.53 


2.37 ± 1.76 


N/N 


A2199 


16: 


:28:37.126 


+39:32:53.29 


0.0302 


2.69 ± 


0.42 


2.64 ± 


0.17 


2.97 


±0.30 


2.26 ± 0.06 


Y/S 


A2255 


17: 


:12;54.538 


+64:03:51.46 


0.0800 


7.13 ± 


1.38 


4.08 ± 


0.26 


7.11 


±0.33 


2.99 ± 1.79 


N/N 


A2589 


23: 


:23;56.772 


+ 16:46:33.19 


0.0416 


3.03 ± 


0.75 


1.88 ± 


0.12 


2.54 


±0.17 


1.77 ±2.04 


Y/W 



Note: The cluster mass, A/sqo.m-m.j, , is derived from the M500 - Mjas.soo relation as tabulated in Lagana et al. (2011), and only used for 
comparison with the cluster mass, M500, derived from the "harmonic" velocity dispersion. "S", "W", and "N" denote strong cool-core, weak 
cool-core, and noncool-core clusters (see Hudson et al. 2010). 

Table 2. Slopes and characteristic magnitudes of the double Schechter function fit to the GLF, gas-mass fraction, and stellar mass 
fraction. 



Name 








M..1 




an 




M,2 






/gas.500 


f. 


..500 


A0085 


-1 


.484 ± 


0.045 


-21.87 ± 


0.71 


-1.945 ± 


0.072 


-20.30 ± 


0.05 


0. 


128 


±0.010 


0.0115 


±0.0012 


A0400 


-1 


.194 ± 


0.057 


-22.24 ± 


0.11 


-1.821 ± 


0.027 


-18.94 ± 


0.06 


0. 


,075 


± 0.008 


0.0240 


± 0.0039 


IIIZw54 


-1 


.289 ± 


0.120 


-22.00 ± 


0.70 


-1.651 ± 


0.028 


-18.42 ± 


0.01 


0. 


076 


±0.013 


0.0239 


± 0.0039 


A1367 


-1 


.184 ± 


0.029 


-22.72 ± 


0.81 


-1.922 ± 


0.062 


-18.98 ± 


0.04 


0. 


118 


± 0.009 


0.0247 


± 0.0028 


MKW4 


-1 


.666 ± 


0.073 


-20.62 ± 


0.35 


-1.802 ± 


0.024 


-19.65 ± 


0.10 


0. 


,094 


±0.013 


0.0233 


± 0.0040 


ZwC11215 


-1 


.698 ± 


0.013 


-22.02 ± 


0.11 










0. 


124 


±0.013 


0.0143 


±0.0017 


A 1650 


-1 


.271 ± 


0.051 


-23.10 ± 


0.60 










0. 


148 


±0.018 


0.0217 


± 0.0026 


Coma 


-1 


.303 ± 


0.025 


-22.34 ± 


0.36 


-1.935 ± 


0.019 


-19.76 ± 


0.08 


0. 


129 


± 0.009 


0.0201 


±0.0018 


A 1795 


-1 


.327 ± 


0.006 


-21.48 ± 


0.61 


-1.847 ± 


0.032 


-20.17 ± 


0.02 


0. 


150 


±0.014 


0.0182 


± 0.0022 


MKW8 


-1 


.145 ± 


0.074 


-20.55 ± 


0.45 


- 1.959 ± 


0.028 


-19.32 ± 


0.02 


0. 


128 


±0.016 


0.0259 


±0.0031 


A2029 


-1 


.227 ± 


0.010 


-21.83 ± 


0.86 










0. 


091 


± 0.008 


0.0065 


± 0.0007 


A2052 


-1 


.313 ± 


0.005 


-21.74 ± 


0.81 


- 1.964 ± 


0.014 


-19.68 ± 


0.01 


0. 


134 


±0.014 


0.0254 


± 0.0030 


MKW3S 


-1 


.422 ± 


0.033 


-22.62 ± 


0.49 


-1.822 ± 


0.036 


-19.81 ± 


0.03 


0. 


147 


±0.017 


0.0268 


± 0.0035 


A2065 


-1 


.061 ± 


0.029 


-22.01 ± 


0.26 


-1.643 ± 


0.015 


-20.46 ± 


0.01 


0. 


069 


± 0.008 


0.0065 


± 0.0006 


A2142 


-1 


.271 ± 


0.014 


-21.87 ± 


0.10 


- 1.495 ± 


0.094 


-20.38 ± 


0.01 


0. 


187 


±0.017 


0.0114 


±0.0011 


A2147 


-1 


.304 ± 


0.009 


-21.52 ± 


0.12 


-1.879 ± 


0.015 


-19.91 ± 


0.01 


0. 


114 


±0.010 


0.0154 


±0.0015 


A2199 


-1 


.154 ± 


0.031 


-21.44 ± 


0.34 


-2.024 ± 


0.014 


-19.65 ± 


0.01 


0. 


110 


±0.010 


0.0177 


±0.0017 


A2255 


-1 


.419 ± 


0.024 


-23.02 ± 


0.98 










0. 


100 


±0.010 


0.0095 


±0.0011 


A2589 


-1 


.085 ± 


0.013 


-21.71 ± 


0.17 


-1.838 ± 


0.021 


-19.89 ± 


0.01 


0. 


,084 


±0.011 


0.0169 


± 0.0023 
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Table 3. Comparison of the best fits of the gas-mass fraction, stellar mass fraction, and baryon-mass fraction as a function of the 
cluster total mass between different samples. 



Sample 


Best fit 




Intrinsic scatter 


19clusters 


/gas,500 ■ 


= 10-'"-''^^*"-^"*(M5oo/[10'"Mo])-»""*""^' 


(26 ± 4)% 


19clusters-A2029-A2065 


/gas,500 ■ 


= 10-"-'«±''i«(M5oo/[10'''Mo]r8±<'-36 


(23 ± 17)% 


19clusters and Sun-l-09 


/gas,500 - 


= 10-"-'"±''»2)(M5oo/[10'''Mo])''-3<'±o.o7 


(26 ± 8)% 


Pratt+09 


^gas,500 ■ 


= 10-"-''2''±»-''""'""W(M5oo/[2 X 10'^Mo])«'-2'*''°3"'""<^' 


(12 ± 2)% 


Sun-l-09 (Tier 1 groups -I- clusters) 


/gas.SOO - 


= (0.0724 ± 0.0078)/7„'-5(M5oo/[10'Vi„'Mo])°'»3±0-03i 




Andreon 2010 


/gas.SOO - 


= 10-<"-'"*''»«(M,oo/[10'*-=Mo])'''5±0-''3 


(14 ± 2)% 


19clusters 


/.,500 = 


10-<'"±0.05)(yi^_^^/[10l4]y[^])-(0.49±0.09) 


(20 ± 4)% 


1 9clusters- A2029- A2065 
Lin-l-03 


/•.500 = 
/•.500 = 


jQ-(1.56±0.05)(^^^^^^0''*Mo])-<°-3'*''°''> 

.nnnm/ r, \ -(0.26+0.09) 

0.0164!°™° (M5oo/[3 X 10»Mo]) 


(15 ± 29)% 


Gonzalez-l-07 


/■.,500 = 


jQ7.57±0.08)y(^^^^-(0.64±0.13) 




Giodini-l-09 


/t,500 = 


(0.050 + 0.001)(M5oo/[5 x 10''Mo])-<''-*±°"'" 


0.35 dex 


Giodini-l-09 and Lin-l-03 


,/*.500 = 


(0.050 + 0.001)(M5oo/[5 x 10"Mo])-<''3''±°'"*' 


0.50 dex 


Andreon 2010 


/*,500 


(M5oo/[10'«Me])-<''«-°-''«' 


(41 ± 5)% 


19clusters 


/b,500 = 


10-'""±»-^'"(M5oo/[10"'Mo])-(»-^''*»-^«» 


(45 ± 6)% 


1 9clusters- A2029- A2065 
Lin-l-03 


/b,500 = 
/b,500 = 


10-(0.96±0.25)(^^^^/[l(3l4jy[^-|-,0.22rf.57 

,„„„-/ , . \0.148±0.040 

0.148:°;°°J(M5oo/[3x lO'^Mo]) 


(29 ± 14)% 


Giodini-l-09 


/b,500 = 


(0.123 ± 0.003)(M5oo/[2 x 10"*Mo])''°''±°''3 
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Fig. 1. Cluster mass derived from the velocity dispersion compared to the cluster mass derived from the gas mass using the mass 
vs. gas mass relation. The solid line shows the best fit. The dashed line denotes the 1:1 ratio. The green circles, red stars, and 
black diamonds stand for noncool-core, weak cool-core and strong cool-core clusters, respectively. The two most outliers are 
A2029 and A2065, which have the highest masses derived from the harmonic velocity dispersion. 
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Fig. 2. Gas-mass fraction (red, open circles), stellar mass fraction (green, open stars), and baryon-mass fraction (black, open 
triangles) as a function of the total mass and the best fits excluding A2065 and A2029 in the same colors with solid line. The 
stellar mass fraction in our observational sample decreases with increasing cluster mass with (20+4)% intrinsic scatter for the 19 
clusters. For comparison, we also show the gas-mass fraction (gray, solid circles), stellar mass fraction in galaxies (gray, solid 
stars), and the sum of these two fractions (gray, solid triangles) for the simulated sample of 21 clusters with AGN feedback in 
Puchwein et al. (2010), as well as the gas mass fractions of the X-ray selected groups in Sun et al. (2009; gray, open boxes) and 
near-infrared selected groups in Dai et al. (2010; gray, open diamonds). The best fit of the baryon-mass fraction as a function of 
the total mass of the observational sample in Lin et al. (2003) is shown in gray with dot-dashed line. The gray band shows the 
Icr measurement from the WMAP 5-year result (Dunkley et al. 2009). A2142 displays a baryon-mass fraction of 0.198 + 0.008, 
which exceeds the WMAP result with a 3cr significance. The best fit of the gas-mass fractions combining our clusters and the 
clusters in Sun et al. (2009) is shown in red with dashed line with (26±8)% intrinsic scatter 
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Fig. 3. Stellar-to-gas mass ratio within rjoo vs. offset between the X-ray flux-weighted center and the BCG position. The symbols 
have the same meaning as those in Fig.[T] The intrinsic scatter is (41+8)% and (66+11)% for the ten high- and nine low-offset 
systems, respectively. This indicates that the scatter of the stellar-to-gas mass ratios may reflect the freedom in the formation 
history of the systems since the last major merging episode. 
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Fig. 4. Stellar-to-gas mass ratio within rsoo vs. gas entropy at 0.2r5oo and the best fit in the upper panel, as well as the residuals 
normalized by the error bars in the lower panel. The symbols have the same meaning as those in Fig. [1] The data support the 
interpretation that heating from merging quenches the star-formation activity of galaxies in massive systems, and feedback from 
supernovae and/or radio galaxies drives a significant amount of gas to the regions beyond r500 or, alternatively, a substantially 
higher stellar mass fraction in the ICL is present in nonrelaxed systems (e.g., Pierini et al. 2008). 
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Fig. 5. Stellar-to-gas mass ratio within rsoo vs. total mass (left panel) and iron abundance (right panel) and their best fits. The 
symbols have the same meaning as those in Fig.[T] These two correlations indicate that the iron in the ICM mainly comes from 
the pollution by the star formation that happened in the past. In less massive galaxy systems, the star formation efficiency is 
higher. 




Fig. 6. Iron abundance (left panel) and iron mass (right panel) vs. total mass and their best fits. The symbols have the same 
meaning as those in Fig. [1] This finding agrees with the trend seen in simulations (e.g.. Figure 1 1 in Fabjan et al. 2010) toward 
less massive clusters having lower gas-mass fractions but higher iron-mass fractions, so that the gas is more metal-rich in those 
systems. 




Fig. 7. Total mass-to-optical light ratio within rgoo vs. total mass and the best fit excluding A2029 and A2065. The symbols have 
the same meaning as those in Fig. [1] The behavior of the total mass-to-optical light ratio supports the variation in star-formation 
efficiencies. 



